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(57) ABSTRACT 

A whole-body mathematical model (10) for simulating 
intracranial pressure dynamics. In one embodiment, model 
(10) includes 17 interacting compartments, of which nine lie 
entirely outside of intracranial vault (14). Compartments (F) 
and (T) are defined to distinguish ventricular from extra- 
ventricular CSF. The vasculature of the intracranial system 
within cranial vault (14) is also subdivided into five com- 
partments (A, C, P, V, and S, respectively) representing the 
intracranial arteries, capillaries, choroid plexus, veins, and 
venous sinus. The body’s extracranial systemic vasculature 
is divided into six compartments (I, J, O, Z, D, and X, 
respectively) representing the arteries, capillaries, and veins 
of the central body and the lower body. Compartments (G) 
and (B) include tissue and the associated interstitial fluid in 
the intracranial and lower regions. Compartment (Y) is a 
composite involving the tissues, organs, and pulmonary 
circulation of the central body and compartment (M) repre- 
sents the external environment. 

35 Claims, 9 Drawing Sheets 
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WHOLE-BODY MATHEMATICAL MODEL 
FOR SIMULATING INTRACRANIAL 
PRESSURE DYNAMICS 

RELATED APPLICATION DATA 5 

This application claims the benefit of U.S. Provisional 
Patent Application No. 60/409,551, filed Sep. 10, 2002, 
entitled “A Mathematical Model for Human Intracranial 
Pressure and Flow Dynamics,” and U.S. Provisional Patent to 
Application No. 60/416,508, filed Oct. 7, 2002, entitled “A 
Mathematical Model for Human Intracranial Pressure and 
Flow Dynamics.” 

GOVERNMENT LICENSE RIGHTS 15 

The U.S. Government may have a paid-up license in this 
invention and the right in limited circumstances to require 
the patent owner to license others on reasonable terms as 
provided for by the terms of Grant No. DMS-96-26391 20 
awarded by the National Science Foundation, and under 
Grant No. NGT5-40110 and Cooperative Agreement No. 
NCC5-581 awarded by NASA. 

FIELD OF THE INVENTION 25 

The present invention relates generally to mathematical 
models for simulating intracranial pressure dynamics and 
more particularly relates to a whole-body mathematical 
model for simulating intracranial pressure dynamics. 30 

BACKGROUND 

Lumped-parameter models represent an attractive method 
for examining pressure dynamics involving complicated 35 
human physiology. In a lumped-parameter modeling 
approach, the physiological system is subdivided into a 
number of linked, interacting subunits termed “compart- 
ments.” In general, each compartment contains a single 
physical constituent, e.g., blood, cerebrospinalfluid (CSF), 40 
or tissue and interstitial fluid. However, depending on the 
model’s complexity, a given constituent may appear in more 
than one compartment of the model. Dynamics in each 
compartment is specified by lumped, time-dependent func- 
tions giving compartmental pressures, while incremental 45 
changes in flows and compartmental volumes are obtained 
by associating resistance and compliance parameters with 
adjacent compartments. In particular, interaction between 
adjacent subunits is assumed to take place at the interfaces 
of the model’s compartments. 50 

With few exceptions, previous models of this type have 
adopted restrictions known as the “Kellie-Monro Doctrine” 
to reduce complexity. The Kellie-Monro framework consid- 
ers the intracranial system to be completely enclosed within 
the intracranial vault, which is assumed to be rigid. A 55 
specified inflow of blood to the intracranial arteries provides 
a forcing for the system, and outflow from the jugular bulb 
is assumed to instantaneously equate to this inflow. These 
restrictions yield a closed system with constant total volume. 
Strictly intracranial models have produced a number of 60 
important results that illuminate the mechanisms of intrac- 
ranial pressure adjustments in situations involving both 
normal and pathophysiology. However, the ability of these 
closed-system models to incorporate the influence of impor- 
tant extracranial factors on intracranial pressure dynamics is 65 
clearly limited. For example, the important buffering effects 
of the spinal CSF space on intracranial pressure cannot be 
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directly included. From a mathematical point of view, the 
constant volume constraint also produces an over deter- 
mined system of equations that requires special handling to 
avoid singular behavior. 

Lumped-parameter compartmental models of the present 
type have a long history, dating to the earliest such model of 
the intracranial system formulated by Monro in 1783. This 
first model was bi-compartmental, considering incompress- 
ible brain matter and blood as its two constituents. In the 
work of Monro’s student Kellie 40 years later, the vascular 
compartment was further subdivided into arterial and venous 
blood to produce a three-compartment model. Since the 
pioneering work of Monroe and Kellie, increasingly more 
complex models of the intracranial compartment have been 
posited in order to more realistically describe the relation- 
ship between intracranial pressures and volumes. There has 
been a steady increase in the number of fluid compartments, 
the introduction of a separate cerebrospinal fluid compart- 
ment, the inclusion of cardiovascular input, and a relaxation 
of the treatment of system constituents as incompressible 
fluids and matter. As noted in prior art studies, the intrac- 
ranial system involves a number of subsystems, which 
interact through complex mechanisms, and the classical 
piecewise approach, which often considers CSF dynamics 
separately, is not suited to studying multiple parameter 
changes and the effects of interconnected subsystems on 
each other. By contrast, lumped-parameter models of the 
intracranial system are capable of including and linking 
different subsystems, so that such interactions can be exam- 
ined. 

When considering lumped parameter models, it is impor- 
tant to realize that a compartment does not necessarily 
correspond to a precise physical location in the body. For 
example, with a single CSF compartment in a model, CSF 
in the ventricles cannot be distinguished from CSF in the 
subarachnoid and spinal spaces. This is one of the main 
limitations of the lumped-parameter approach. Additional 
spatial resolution can be realized only by subdividing the 
physical system into a larger number of compartments based 
on spatial considerations. For example, distinct ventricular 
CSF and extra-ventricular CSF compartments may be 
included as opposed to a single lumped CSF compartment. 
In principle, the entire body could be finely subdivided in 
this manner into separate compartments to provide the 
desired degree of spatial resolution. However, clearly this 
subdivision process cannot be carried to an extreme as the 
resulting system of linked governing equations will rapidly 
become too large for practical analysis and solution. 

Despite their evolving complexity, two common features 
characterize most earlier lumped parameter models for pres- 
sure dynamics in the intracranial system. The first common 
feature is an assumption that all resistance and compliance 
parameters can be represented by constants. This leads to a 
linear system of governing differential equations. The sec- 
ond common feature is adoption of the “Kellie-Monro 
Doctrine,” which assumes that the intracranial system can be 
confined within the cranial vault. By requiring that inflow to 
the intracranial arteries equals outflow from the jugular bulb, 
this assumption produces a closed system that conserves 
total intracranial volume. 

As mentioned above, when the intracranial space is 
treated as a closed volume-conserving system contained 
within the (nearly) rigid cranial vault, important mecha- 
nisms for the influence of extracranial physiology on intrac- 
ranial pressure dynamics cannot be included in the resulting 
models. For example, the ability of the spinal portion of CSF 
space to buffer fluctuations of intracranial CSF pressures 
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cannot be directly introduced under the Kellie-Monro Doc- 
trine. At least two prior art mathematical models of intrac- 
ranial pressure dynamics include aspects of extracranial 
physiology. The first model includes a compliance for CSF 
storage within the lumbar channel. The first model contains 
three compliances, four resistances, and involves differential 
equations based on a hydrodynamic model for the physical 
system and its electrical circuit equivalent. The first model 
allows the dynamic relationship between cerebral perfusion 
pressure, intracranial pressure, and cerebral blood flow in 
various states of autoregulation to be studied. Use of the first 
model in conjunction with clinical data has determined 
which indices can be derived using transcranial Doppler 
ultrasonography and which trends of intracranial pressure 
and blood pressure are useful in clinical tests of autoregu- 
latory reserve. However, despite not strictly abiding by the 
Kellie-Monro Doctrine, the first model falls far short of 
being a whole-body model for intracranial pressure dynam- 
ics. While the first model can include direct buffering effects 
of CSF within the lumbar channel, it does not contain other 
important extracranial physiology. For example, the first 
model fails to include representations for the extracranial 
physiology of the cardiovascular system and the tissues and 
extracellular fluids of both the central and lower body, 
including cardiovascular autoregulation, colloid osmotic 
pressure, and a lymphatic system. 

The second model uses a single ground compartment to 
represent the portion of the body below the clavicles. The 
second model contains three resistances, including a resis- 
tance between the intracranial veins and the rest-of-body 
compartments. Thus, outflow from the intracranial system 
depends on pressure differences and need not instanta- 
neously equal the specified inflow. The second model also 
contains three compliances, including a compliance between 
the CSF and rest-of-body compartment that represents the 
ability of the spinal portion of CSF space to expand and 
buffer CSF pressures by compressing the large veins that 
surround the spinal theca. Two of the three model compli- 
ances are pressure dependent. The associated functional 
forms for the pressure-dependent compliances determined in 
the second model are used in the present work. The second 
model still falls short of being a whole-body model for 
intracranial pressure dynamics. The second model merely 
lumps extracranial physiology into a single compartment, 
rather than finely subdividing it into a plurality of compart- 
ments. In addition, the second model fails to include the 
external environment as an implicit compartment. 

SUMMARY OF THE INVENTION 

One aspect of the present invention is a method of 
modeling pressure dynamics of a body’s intracranial system 
including the following steps: dividing the body into a 
plurality of compartments and a heart pump, each of the 
plurality of compartments representing a portion of the body, 
the heart pump interacting with at least one of the plurality 
of compartments; deriving a plurality of differential flow 
equations, each of the plurality of differential flow equations 
corresponding to one of the plurality of compartments; and 
solving the plurality of differential flow equations. 

Another aspect of the present invention is a system for 
modeling an intracranial system including a body compart- 
ment module adapted to divide the body into a plurality of 
compartments and a heart pump, each of the plurality of 
compartments representing a portion of the body, a flow 
equation module adapted to derive a plurality of differential 
flow equations, each of the plurality of differential flow 
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equations corresponding to one of the plurality of compart- 
ments, and an equation solver module adapted to solve the 
plurality of differential flow equations. 

Still another aspect of the present invention is a method of 
5 modeling pressure dynamics of an intracranial system that 
further accounts for cerebrovascular autoregulation by a 
sympathetic nervous system. 

Yet another aspect of the present invention is a method of 
modeling pressure dynamics of a body’s intracranial system 
to including the following steps: dividing the body into a 
plurality of compartments and a heart pump, each of the 
plurality of compartments representing a portion of the body, 
the heart pump interacting with at least one of the plurality 
of compartments, wherein a plurality of the plurality of 
15 compartments are vascular and a plurality of the plurality of 
compartments are non-vascular, the vascular compartments 
including at least one of the intracranial arteries, intracranial 
capillaries, choroids plexus capillaries, venous sinus jugular 
veins, intracranial veins, central arteries, central capillaries, 
20 central veins, extra-ventricular CSF, lower arteries, lower 
capillaries, and lower veins, the non-vascular compartments 
including at least one of lower tissue, brain, ventricular CSF, 
and extra-ventricular CSF ; deriving a plurality of differential 
flow equations, each of the plurality of differential flow 
25 equations corresponding to one of the plurality of compart- 
ments; and solving the plurality of differential flow equa- 
tions. 

BRIEF DESCRIPTION OF THE DRAWINGS 

30 

For the purpose of illustrating the invention, the drawings 
show a form of the invention that is presently preferred. 
However, it should be understood that the present invention 
is not limited to the precise arrangements and instrumen- 
35 talities shown in the drawings, wherein: 

FIG. 1 is a simplified schematic of one embodiment of a 
multi-compartmental model for determining intracranial 
pressure dynamics; 

FIG. 2 is a block diagram of the various compartments of 
40 a model according to FIG. 1; 

FIG. 3 is the intracranial arteries compartment of the 
model according to FIG. 1; 

FIG. 4 is a chart of interstitial fluid pressure vs. relative 
lymph flow; 

45 FIG. 5 is a chart of percent reduction in Q AC VS. P^; 

FIG. 6 is a simplified schematic of a computer system 
including one embodiment of the present invention; 

FIG. 7 A is a chart of pulsatile cardiac forcing (Q^) over 
time; 

FIG. 7B is a chart of the 120/80 central artery pressure 
response over time (solid) and the 100/65 intracranial artery 
pressure response (dashed); 

FIG. 8 A is a chart of the predicted CSF pressure vs. 
55 volume; 

FIG. 8B is a chart of conductance to CSF output accord- 
ing to one embodiment of the present invention; 

FIG. 9 is a chart of the percentage of Q^ c drop due to a 
positional change with all SNS reflexes activated (solid) and 
60 without any SNS reflexes included (dashed); 

FIG. 10 is a chart of percentage decrease of P 7 (solid) and 
of Q^ 7 (dashed) as a function of percentage blood loss due 
to hemorrhage; and 

FIG. 11 is a chart of the predicted percentage of Q^ c drop 
65 vs. percentage of blood loss due to hemorrhage with both 
SNSz reflexes activated (solid) and no SNSz reflexes 
included (dashed). 
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DETAILED DESCRIPTION OF THE 
INVENTION 

The present invention is directed to a whole-body math- 
ematical model for simulating intracranial pressure dynam- 
ics. Examples of the disclosed invention are depicted in 
FIGS. 1-11, although it should be understood that the 
present invention is not limited to this (or any other) 
particular embodiment, but rather is intended to cover mod- 
els and systems that fairly fall within the broad scope of the 
appended claims. 

As described inU.S. Provisional Patent Serial No. 60/409, 
551, which is incorporated by reference as if disclosed 
herein in its entirety, the model of the present invention 
revokes the Kellie-Monro Doctrine by consistently embed- 
ding the intracranial system within whole-body physiology. 
In one embodiment, the model of the present invention 
includes 17 interacting compartments, of which nine lie 
entirely outside of the intracranial vault. Two distinct com- 
partments are defined to distinguish ventricular from extra- 
ventricular CSF. The vasculature of the intracranial system 
within the cranial vault is also subdivided into five com- 
partments representing fluid in the intracranial arteries, 
capillaries, choroid plexus, veins, and venous sinus. The 
body’s extracranial systemic vasculature is divided into six 
compartments representing the arteries, capillaries, and 
veins of the central body and the lower body. Tissue and the 
associated interstitial fluid in the intracranial and lower 
regions are divided into two compartments. A composite 
compartment involving the tissues, organs, and pulmonary 
circulation of the central body and an implicit compartment 
representing the external environment complete the model. 
Since the time-dependent compartmental pressure functions 
are obtained from physical pressures through a “lumping” 
procedure that involves space-averaging over the physical 
extent of a compartment, the subdivision of physical con- 
stituents into distinct spatially limited compartments is nec- 
essary to provide spatial resolution in this modeling 
approach. 

In addition to allowing direct flows (e.g. arteries to 
capillaries to veins) between adjacent compartments, the 
present work includes the extracranial transfer of fluid 
between capillaries and tissue by filtration. An extracranial 
lymphatic system is also included in the model. Components 
of the model allow regulation of systemic vascular pressures 
by the sympathetic nervous system, and, at less than extreme 
(high or low) pressures, autoregulation mechanisms provide 
constant blood flow to the cerebrovascular capillaries and 
the choroid plexus as well as constant production of CSF in 
the choroid plexus. Fluid intake, renal output of fluid, and 
adjustment of body volume in response to changes in 
ambient environmental pressure are allowed. A realistic 
representation for cardiac uptake and cardiac output pro- 
vides the forcing for this system. 

The form of the present model is a system of governing 
differential equations for the fully time-dependent compart- 
mental pressure functions. A component of the present work 
is determining appropriate forms for the non-constant resis- 
tance and compliance parameters in the model, which may 
be functions of both pressures and time. Calibration of 
physically realistic scale values for parameters and flows is 
also a step in the simulation process. 

Turning to the drawings, wherein like reference numerals 
refer to like elements, FIG. 1 illustrates one embodiment of 
a mathematical model 10 according to the present invention. 
In mathematical model 10, the human body is divided into 
16 distinct compartments, lettered A, B, C, D, F, G, I, J, O, 
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P, S, T, V, X, Y, and Z. Compartment A represents the 
intracranial arteries. Compartment B represents the brain. 
Compartment C represents the intracranial capillaries. Com- 
partment D represents the lower capillaries. Compartment F 
5 represents the ventricular CSF. Compartment G represents 
the lower tissue. Compartment I represents the central 
arteries. Compartment J represents the central capillaries. 
Compartment O represents the central veins. Compartment 
P represents the choroid plexus capillaries. Compartment S 
to represents the venous sinus jugular veins. Compartment T 
represents the extra- ventricular CSF. Compartment V rep- 
resents the intracranial veins. Compartment X represents the 
lower veins. Compartment Y represents the rest of the body. 
Compartment Z represents the lower arteries. Compartment 
15 M represents the external atmosphere. Compartment H L and 
Hr represent the left and right heart chambers. The physical 
constituents in subunits of the present model are blood, CSF, 
and tissue and interstitial fluid. With a single exception, i.e., 
compartment Y, each compartment is composed of a single 
20 constituent. Spatial resolution is obtained by first dividing 
the body into intracranial and extracranial components. To 
help delineate the intracranial system, a thick line 12 in FIG. 
1 represents a cranial vault 14, which differentiates intrac- 
ranial from extra-cranial compartments. Compartments I, Z, 
25 Y, J, D, G, O, and X lie completely outside of cranial vault 
14, and compartment T includes both intracranial and 
extracranial physiology. 

As illustrated in FIG. 2, the majority of the compartments 
in the current model are vascular. The 11 vascular compart - 
30 ments can be subdivided spatially into three groups: intrac- 
ranial compartments A, C, P, V, and S; central compartments 
I, J, and O, and lower compartments Z, D, and X. The 
“lower” compartments represent the region below the pelvis. 
The “central” compartments in the model represent the 
35 region between the lower body and the clavicles and also 
include extracranial body components above the clavicles 
but outside of the cranial wall. The vascular system in each 
of these three regions, i.e., intracranial, central, and lower, is 
divided into separate artery, capillary, and vein compart- 
40 ments. However, in the intracranial space, the choroid 
plexus capillaries are placed in a separate compartment P 
from the rest of the intracranial capillary bed. This allows 
these capillaries to autoregulate to maintain the production 
of CSF when ventricular CSF pressures rise or fall. The 
45 venus- sinus veins where CSF absorption occurs through the 
arachnoid villa and the jugular veins are also placed in a 
compartment separate from the remainder of the intracranial 
veins, i.e., compartment S. There are four strictly non- 
vascular model subunits. Two of these compartments rep- 
50 resent tissue matter and interstitial fluid, i.e., lower com- 
partment G and brain compartment B, and the other two 
represent CSF compartments, i.e., ventricular F and extra- 
ventricular T. Compartment T is both an intracranial and a 
central compartment. Compartment T, which contains extra- 
55 ventricular CSF and lies partially outside of the cranial vault, 
includes both the subarachnoid and spinal portions of CSF 
space. Compartment T serves as a bridging compartment 
and explicitly allows for buffering of CSF pressures by the 
spinal theca. Two of the three regions contain a separate 
60 compartment for the tissue and interstitial fluid. The excep- 
tion is the central region where the tissue matter, organs 
(except for the heart), interstitial fluid and pulmonary cir- 
culation are lumped into a composite rest-of-body compart- 
ment Y. The central region contains an explicit representa- 
65 tion for the heart pump having a left chamber pump H z , and 
a right chamber pump Hr, and a realistic pulsatile cardiac 
output provides the major forcing in the current model. 
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External environment compartment M is also an implicit 
17-th subunit in the model. No attempt has been made to 
depict relative volumes in FIG. 1, and hence the relative 
sizes of these compartments in the figure do not reflect 
relative volumes. 

The pressure dynamics of the intracranial system are 
governed by a system of differential equations within math- 
ematical model 10. Four basic assumptions lead to these 
equations: 

(1) all fluids are considered incompressible and isother- 
mal; 

(2) pressure driven flows are laminar and related to 
pressure differences by 


Qu = 



= Z u {P i -P j ) = Z u P u , 


(l) 


where is the flow from compartment i into compartment 
j, P z and Pj are the spatially -averaged pressures of compart- 
ments i and j respectively, R^. is the lumped resistance, Z tj is 
the fluidity (inverse of R^), the pressure difference P^P^-P^, 2 5 
and R^=-R^; 

(3) in the case of fluid filtration from the capillaries into 
the interstitial space, the flow is governed by the Starling- 
Landis Equation, i.e., 

30 

Yi\Mio^=K c l{P-P t )-o c l% c -%^)=K c lP c -o c pi ct \ ( 2 ) 

where P c is the capillary pressure, P t is the interstitial fluid 
pressure, jc c is the blood plasma colloid osmotic pressure, n t 
is the interstitial fluid colloid osmotic pressure, K ct is the 35 
filtration coefficient, and a ct is the capillary membrane 
reflection coefficient. The notation for pressure difference 
has been extended to osmotic pressure differences by defin- 
ing 7t ct =7Z c -7i t ; and 

(4) the deformation of the membrane between adjacent 40 
compartments is a function of the change in pressure dif- 
ference between these compartments, i.e., 


dVy d(Pi-Pj) ^ d(Pij) (3) 

dt ~ ij dt ~ ij dt 


where V,, denotes the instantaneous volume of the ‘cup’ 
formed in the membrane at the interface of compartments i 
and j, C^- denotes the compliance between these two com- 
partments, and C^C^. 

As follows, the system of governing equations is obtained 
by invoking the conservation law 55 

flow rate in-flow rate out=rate of volume change (4) 


8 

Specifically, the pressure dynamics of compartment I are 
governed by the following equation: 


Qhi ~ {Z,jP,j + Z/zPiz + Z m P m ) = C it 


d P/t 
dt 


dVy/ 

~dt~' 


(5) 


Note that in equation (5), the subscript “IY” has been 
reversed thereby causing the associated term to be negative. 
Also, in both FIGS. 1 and 3, a filled arrow indicates a 
one-way flow and a hollow arrow indicates a pressure 
dependent resistance. In terms of pressure differences, the 
relation of equation (4) yields the following additional 12 
differential equations and three scalar equations for the 
remaining 15 body compartments: 

Central Capillary Compartment (J): 

Z/jP ij- (K jy (P JY - VjyXjy) +Z JO P jo) =0 ; (6) 

Rest of Body Compartment (Y): 


Qw + Kjy(Pjy - & jy^jy) - ( Qyo + Qk ) = 


(7) 


dP Y M dV Y l dVyo 

Cm ~dT + ~dr + ~dT' 


Central Venous Compartment (O): 


ZjqP jo + ZxqPxo + Z$qP so + Zro P to + Qco + Qyo - Qoh = 


( 8 ) 


dP 0 T 
7 dt 


dV Y p _ 

dt 


Lower Artery Compartment (Z): 


Z/zP/z - ZzdPzd = Czc - 


dPzc 

dt 


(9) 


Lower Capillary Compartment (D): 

ZzdP zd~ (Zdx? dx*Kdg(P dg~ g dg ji dg)) = ^ : > 

Lower Tissue Compartment (G): 


( 10 ) 


Kdg(Pdg - cr dg^dg) - Qgo = 


(11) 


„ dP G z „ dP G x „ dPcM 
Czc ~dT + Cgx ~dT + Cgm ~dt~’ 


Lower Venous Compartment (X): 


ZdxPdx -ZxqPxo = Cgx 


dPxc 

dt 


(12) 


in each compartment. Referring now to FIG. 3, compartment 
I from FIG. 1 is illustrated. Compartment I represents the 
central arteries. The flow rate into compartment I is repre- 
sented by Qjtj from left heart chamber pump H L . The flow 
rate out of compartment I is the sum of the flows associated 


Intracranial Artery Compartment (A): 


60 


Z/aPia ~ ( ZacPac +Z A pP ap) = C A b 


dP ab 
~dt ' 


• C A f 


dP af _ 
~dt 


(13) 


equation (1). The rate of volume change of compartment I is 
equal to the sum of the terms associated with the compli- 
ances C IT and C IY , as represented by equation (3). 


65 


Intracranial Capillary Compartment (C): 

ZacP act(Z cb P cb +Z cv P cf) = 0; 


(14) 
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Choroid Plexus Compartment (P): 


Zap Pap ~ (Z PF P PF + Z PV P PV ) = C PF 


dP PF 
dt ’ 


Intracranial Veins Compartment (V): 


ZcvPcv +Z PV P PV +Z BV P BV -Z vs Pvs = 


_ dPvB dPvF 
CBV — +CFV —’ 


( 15 ) 


10 

The set of governing equations may be summed to derive 
the following constraint regarding compliances between the 
body and the external environment: 


d Pym d Pgm 

(Qw ~ Qk) + ( Qhi ~ Qoh ) = Cym 


(21) 


If cardiac output equals cardiac uptake (Q^ 7 =Q 0i7 ) and 
the atmospheric pressure is constant, equation (21) simpli- 
fies to 


Venous Sinus — Jugular Veins Compartment (S): 


15 


dP Y dP G 
Qw - Qk = Cym -jj- + Cgm ■ 


(22) 


ZvsPvs + ZtsPts - ZsoPso - Cts 


dPsr 
dt ’ 


( 17 ) 


Ventricular CSF Compartment (F): 


Zp F Ppf - (Z FB P FB + ZptPft ) = 


( 18 ) 


d P fa dP FP dP FB dP F y 

c " -nr +c "^r +c ™— +c "— ; 


Extra- Ventricular CSF Compartment (T): 


Z F T P FT + Z B J P B J - (Zfs PtS + Zjo Pro ) = (19) 

dPjs dPjB dPjo dPn 
C„— + C B J— + C T0 — + C IT —-, and 


Brain Compartment (B): 


If, in addition, Q^Q^M) or fluid intake equals fluid 
discharge (Q^Q^), then equation (22) implies as expected 
20 that the net volume change of the entire body system must 
be zero. Furthermore, if Q^Q^, as will initially be the case 
when a glass of water is consumed, the net volume of the 
body will increase. By equation (22), one of the internal 
compartmental pressures: P G , or most likely P r , must now 
25 increase. However, if as will be the case, C w is large, the 
resultant pressure increase in compartment Y due to the 
volume input from the ambient environment will be small. 

The fluidities in equation (1) must be derived, the filtra- 
tion and reflection coefficients in equation (2) must be 
30 calibrated, and the compliances in equation (3) must be 
calibrated before model 10’s governing differential equa- 
tions may be used in simulations. 

Each constant fluidity Z tj between arbitrary adjacent com- 
35 partments i and j may be calculated by substituting the mean 
pressures and the mean flow rates into equation (1) and 
solving for Z tJ , giving 


ZcbPcB + Zfb P FB - ( Z B yP B y + Z BF P B T ) = 


( 20 ) 40 


Cab- 


dP B A 

dt 


- C B y - 


dP bv 
dt 


- C FB - 


dP BF 

dt 


b C BF - 


dP B j 
dt 


Zy = 


_gtf_ 

Pi-Pj' 


( 23 ) 


The terms dV ^dt and dV ro /dt in the conservation equa- 45 
tions for compartments Y, I, and O have been left intact in 
equations (5), (7), and (8). These volume changes include 
components that reflect the regulation of arterial pressure by 
the sympathetic nervous system. They are considered further 
below. The terms involving dP^/dt in the conservation 50 
equations for 6 compartments Y, and G (equations (7) and 
(11)) denote a volume change for which there is no com- 
pensation through an equal but opposite volume change in 
an adjacent body compartment. Specifically, these terms 
reflect volume increases (or decreases) into the ambient 55 
environment, which is considered infinitely large and unaf- 
fected by pressure changes within the body. Ambient envi- 
ronment compartment M is similar to the ground in an 
electrical circuit model. It will be maintained at a constant 
pressure value in exactly the same way that the ground 60 
voltage in a circuit model is given a prescribed value relative 
to which all other values are gauged. Accordingly, the 
ambient environmental pressure is fixed here at a con- 
stant reference value of zero mmHg. If this is not appropriate 
for a situation being studied, or if the ambient pressure 65 
changes with time, then must be considered as a forcing 
term in the system of equations. 


Therefore, once the mean flows and pressures are esti- 
mated from physical data, the associated constant fluidities 
may be calculated from equation (23). Some model fluidities 
are pressure dependent. A discussion of appropriate expres- 
sions for pressure dependent is provided below where CSF, 
cerebrovascular, sympathetic nervous system, and cardiac 
autoregulation mechanisms are modeled. However, even for 
a non-constant fluidity, a mean scale value may still be 
calculated from equation (23). In the calibrations that follow, 
mean flows and pressures reflect physiological values for an 
average human in a lying down (supine) position. 

To estimate mean flows, it will be assumed that compart- 
mental volumes remain constant in the mean state. Com- 
partmental pressures are certainly pulsatile, so the system 
will have a mean state, but not a steady state. However, since 
volume changes are related to pressure differences between 
adjacent compartments, if pressures in adjacent compart- 
ments move in tandem in the mean state, volumes will 
remain unchanged. Therefore, for a given compartment, if 
all but one of the mean flows are known, the final mean flow 
may be determined from maintaining constant volume. 
Further, once mean flows into and out of a given compart- 
ment are determined, these values provide data for flows in 
the adjacent compartments. While many flows can be esti- 
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mated from data available in published literature, most mean 
flow calibrations must make use of the constant volume 
assumption. 

As a starting point for the calibration of mean flows, 
consider the percentages of cardiac output Q ^ 7 that exit the 
central arteries into the three peripheral artery systems: A, I, 
or Z. These percentages, in decimal form, are given below. 


p^=0.15,=% of cardiac output into intracranial 

region (24) 

p^O.25, 0.35,=% of cardiac output into lower 

region (25) 

Pij=l -(Pi a +Piz) ,=% °f cardiac output into rest of 

body (26) 

Additional helpful percentages and ratios include 

p^0.70,=% of CSF formation from the choroid 

plexus (27) 

X=250,S ac /Q ap (28) 

p IS =0.80,=% of CSF drained into venous sinus (29) 


Literature values give the following mean flows in 
ml/min: 

0^7=5000, 6600, 6900, Cardiac output (30) 

Q cf = 2, Total capillary filtration (with 2 A from liver 

and intestines) (31) 

_ (2 1 pu ) . . (32) 

Qjy = Q cf \ — + — I, Total capillary filtration 

V 3 3 pij 4- pi z ) 

into central tissue 


QdG - O “ ~~QcF-> Total capillary filtration 
3 Pij + Piz 

into lower tissue 


Lower Body Mean Flows 

Qzd = Q/z fr° m constant volume in Z (43) 

Qoo=Q»g f rom constant volume in G (44) 

~Qdx^Qzd~Qdg from constant volume in D (45) 

Qxo = Qdx f rom constant volume in X (46) 

io Intracranial Mean Flows 


Qap = Qia!Q^+ 1) f rom constant volume in A and 

equation (28) (47) 

Qac = ^Qap f rom equation (28) (48) 

15 

Qp^P^/Qf from equation (27) (49) 

Qcb = 0--Pp/)Qf+Qbv f rom Qcb~ Qb r=( 1 ~P P f) Q f (50) 

^ Qcv^QactQcb from constant volume in C (51) 

Qpv = Qap~Qpf from constant volume in P (52) 

Qvs = Qcv + Qpv + Qbv from constant volume in V (53) 

25 Qbt = Qcb + Qfb~Qbv from constant volume in B (54) 

Qft~Qpf^Qfb from constant volume in F (55) 

~Qts = P tsQ,F fr° m equation (29) (56) 

30 

~Qto^~P ts) Qf from constant volume in T (57) 

~Qso from constant volume in S (58) 

Equation (49) states that 70% of the CSF formation comes 
35 from the choroid plexus and equation (50) states that the 
remaining CSF formation comes as a filtrate from the 10 
capillaries through the brain via the Virchow -Robins system 
minus the amount reabsorbed at the venule side of the 
capillaries (Q^). Notice that, as should be the case, 

40 _ _ _ 

Qbt+Qfi—Qf- 


Qf^O. 35, Total CSF formation (34) 

Qsj/= 0.001, Imbalance of diffusion on venule side of 

cerebrovasculature (35) 

0^3=0.044, Flow of CSF through the Virkow-Rob- 

ins Space (36) 

All of the remaining mean flows in the model may now be 
calibrated based on equations (24) through (36) by invoking 
the constant volume assumption in each compartment during 
the mean state. The calibration procedure yields the mean 


flows: 

Central Body Mean Flows 

QoiFQhi from cardiac output=cardiac input (37) 

QiA^iaQifi from equation (24) (38) 

Qir=PizQm from equation (25) (39) 

Qu= Qht~ ( Qia + Qiz) from constant volume in / (40) 

Qjo = Qij~Qjy from constant volume in J (41) 

Qyo=Qyy from constant volume in Y (42) 


In addition to allowing direct flows (e.g. arteries to 
capillaries to veins) between adjacent compartments, the 
present model includes the transfer of fluid between capil- 
45 laries and tissue by filtration. These flows are governed by 
the Starling-Landis equation and are driven by differences 
between the colloid osmotic pressures of the blood plasma 
in the capillaries and the interstitial fluid as well as by the 
usual compartmental pressure differences between capillar- 
50 ies and tissue. Filtration mechanisms are included between 
the capillaries and tissue compartments in the central and 
lower portions of the body in the present model. In the 
intracranial region, significant colloid osmotic pressure dif- 
ferences do occur between the intracranial capillary and 
55 tissues. However, the endothelial cells that make up the 
intracranial capillary wall are so tightly joined that not even 
water molecules can usually pass between them. Thus, under 
normal conditions, colloid osmotic pressure forces in the 
intracranial region are irrelevant. In the case of highly 
60 elevated capillary pressure it is possible for the intracranial 
capillary wall to expand enough to allow water molecules to 
pass between the endothelial cells, and at this point both 
colloid osmotic pressures and regular pressures start gov- 
erning filtration rates. At these elevated pressures, the 
65 description of intracranial capillary filtration as a function of 
capillary pressure will clearly be nonlinear. However, to 
simplify the equation for intracranial filtration, this relation- 
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ship may still be linearly approximated with the slope 
defined by Z CB , Z BV , and Z FB . 

There are two locations in the model where fluid flow is 
governed by the Starling-Landis equation (2). These flows in 
the central and lower body are 

Qjy=K jy ((Pj-P r)-o>y(jtj-Jty)) and (59) 

Qdg = Kdg((P d~P g)~ g dg&d~ :ji g))> (60) 

denoting the flow from the capillaries into the tissue subunits 
of the central and lower body regions, respectively. Solving 
these relations for the filtration coefficients in the mean state 
implies 


i(Pj -Py)-o-jy(xj - Try)) 


((Pd ~ Pc) ~ c dg(xd ~ n c )) 


In this section, the various compliances in equation (3) 
that relate volume adjustments to pressure differences will 
be calculated. In the intracranial region, compartmental 
volume increases are restricted by the rigid cranial wall. 
5 Consequently, compartmental compliances must be pressure 
dependent and diminish to zero as pressure differences 
become large. On the other hand, in extracranial regions, to 
lowest order compliances may consistently be considered 
constant, approximating a linear relationship between pres- 
to sure differences and volume adjustments. The present intrac- 
ranial pressure-dependent compliances are extensions of 
those derived in the prior art, while the constant extra-cranial 
compliances will be derived from estimations of the volume 
and distensibility of each compartment. 

15 In a simplified four-compartment whole-body model for 
CSF dynamics described in the prior art, there are only two 
pressure difference dependent compliances. They allow vol- 
ume adjustments between the CSF and arterial blood and 
between the CSF and venous blood. These compliances have 
20 the general form described by the relation 


The mean flows on the right hand side of these equations 
have been calculated above. Thus, to determine the filtration 


scale values for the mean colloid osmotic pressures and the 
reflection coefficients. 

Mean interstitial fluid colloid osmotic pressure is given by 

G =8 mtnHg (63) 


where P^P^-P^, and the subscripts i and j take the values 
C a f (for arterial/CSF compliance) and C^ 4 (for CSF/venous 
compliance). For both pairs of index values, Coefficients and 
parameters in equation (67) are given by 


and blood plasma colloid osmotic pressure by 

Jt y =Jt z) = 28 mmHg. (64) 

Notice that these values are invariant with respect to 35 
central or lower regions. This is not the case, however, with 
the reflection coefficients o JY and o DG . The reflection coef- 
ficient 11 in the legs is estimated to be approximately 0.9 
while the coefficient of the upper body is less than this value. 
This is reflected by the assignments 40 


The filtration coefficients K JY and in equations (61) and 

(62) are now readily determined. 

Calibration of scale values for all resistance and compli- 
ance parameters from available physical data and other 
relationships must be accomplished before the model’s 50 
governing differential equations can be used in simulations. 

A key step in calibrating model compliances is determining 
the distensibilities of the extracranial compartments of the 
model. In particular, each extracranial compartment has an 
associated total compliance equal to the product of its 55 
associated volume and distensibility. It should be noted that 
calculating compliances by this technique yields a total 
central artery compliance C 7 =1.529 ml/mmHg, which is 
within 5 percent of the data value of 1 .445 ml/mmHg for 
total arterial compliance measured in the ascending aorta by 60 
Chemla et al. The central venous compliance calculated by 
these methods yields C K =50 ml/mmHg while known sys- 
temic venous pressure volume curves suggests a total 
venous compliance of 50 ml/mmHg. Finally, Noordergraaf 
states that the systemic arterial compliance is between 1 and 65 
2 ml/mmHg and systemic venous is between 50 and 200 
ml/mmHg. 


This four-compartment model contains an additional con- 
stant compliance between the CSF and the rest of the body 
(g), which is approximated as 

C^ 4 =0. 13333. (70) 

This feature represents the interface of extra-cranial CSF 
in the spinal theca with the rest of the body. It also acts as 
a background compliance so that the total CSF compliance 
can never reach machine zero in simulations. 

In the present 16 compartment model, the division of the 
body’s cerebrospinalfluid (CSF) space is considerably more 
refined. Hence, the three CSF compliances in the simple four 
compartment model of the prior art must be appropriately 
apportioned among the compliances of the present model. 
Three decimal percentages may be introduced to describe 
this allocation of C^ 4 : 

p^,32 0.164 percentage of C^ 4 allocated to C FV Cl) 

p to =0.214 percentage of C fv 4 allocated to C TO (72) 

p te =0.622 percentage of C^ 4 allocated to C TS (73) 

Equations (71-73) reflect the fact that total CSF volume 
is approximately 140 ml, 23 of which is found in the 
ventricles, 30 in the spinal cord subarachnoid space (theca) 
and the remainder in the cerebral cisterns and subarachnoid 
space. Thus, if the distensibility of the membrane is similar 
in these three different components of CSF, then c > is 
23/140=0.164 of C /v 4 as indicated by p /v in equation (71). 
This same technique yields the values in equations (72) and 
(73). The compliance C PF is excluded in this division of C /v 4 
as the choroid plexus capillaries are known to dilate and 
constrict in order to maintain a constant pressure difference 
between these capillaries and ventricular CSF. This main- 
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tains a constant generation rate of CSF from the choroid 
plexus at normal pressure values. Consequently, a value for 
the compliance C PF is irrelevant in the absence of a pressure 
difference change. 

The ratios that lead to equations (71)-(73) also imply the 
following percentages describing the allocation of C a 4 : 

p a /=0JS6 percentage of C a f allocated to C AF (74) 

p zz =0.214 percentage of C a 4 allocated to C IT (75) 

The background compliance Cf s4 is now divided between 
venous and arterial interfaces based on systemic venous 
volume being four times that of arterial volume, giving 

C v _0.8 C* 4 (76) 

Ca^iafO.2 c fg 4 (77) 

The CSF-related compliances may now be calculated 
based on the above percentages as: 

ml/mmHg 

C ra (P ra )<l.05p A (C/(P ra ) + C^)<: ra =0.036255 

ml/mmHg 

Cts(Pts)= 0.95 p ^1.27626 

ml/mmHg 

C BI (P BT )=0.05 P a (C/(P BT )+C veno J C flI KU37057 
ml/mmHg 

C T o(PTo)=P t o(CA P ro)+C venous )-C TO = 0.200936 
ml/mmHg 

0.0261999 

ml/mmHg 
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and may be experimentally determined by measuring the 
pressure change that is induced by an estimated volume 
change and then taking the inverse of this ratio. Unfortu- 
nately, clinical data does not exist for determination of all of 
5 the compliances required for the current model. However, 
there is information in the literature regarding volumes and 
distensibilities. These are related to compliance by 

Compliance=Distensibility- Volume. (7 9) 

10 Therefore, total compartmental compliance may be cal- 
culated from total compartmental volume (denoted V.) and 
compartmental distensibility (denoted D z ) through equation 
(79). However, before calculating compartmental volumes 
and distensibilities it is useful to describe how the resulting 
15 total compartmental compliances will be allocated to the 
local inter-compartmental compliances. 

Estimations of the local compliances are made by requir- 
ing that the sum of the relevant local compliances must equal 
the total compliance being apportioned. With this restriction, 
20 the model’s predicted compartmental pressure changes 
induced by known volume changes should agree with 
experimental data. Based on this principle, the apportion- 
ment of total compliances among local compliances is 
described by 

25 Central, Local Compliances: 


C I f=M2^[C i -'C it ,0] ml/mmHg 

(80) 

C ro =Max [C 0 - (T 70,0] ml/mmHg 

(81) 

30 C KW =Max[C y -(C /r fC ro ),0] ml/mmHg 

(82) 

Lower, Local Compliances: 


Czg = C-z ml/mmHg 

(83) 


QiiPjr^CjiPjJ+C^J C lf = 0.00571427 
ml/mmHg 

where again, the pressure difference P z -P^ is denoted P^ and 
the scale value 'C i j==C i j(P i -Vjy=C i j(P i j). The compliance C BT , 
introduced to incorporate volume adjustments between the 
brain and subarachnoid CSF, is allocated 5% of C TS . This 
states that of the bulk intracranial compliance in compart- 
ment T, 95% is allocated to the interface with the highly 
distensible venous sinus veins and only 5% to the interface 
with the less distensible brain tissue. A similar allocation is 
made for the bulk intracranial compliance of the ventricular 
CSF compartment F with respect to the intracranial veins 
and the brain tissue. 

Finally, brain/vascular compliances C^(P^) and C BV 
(P BV ) are defined by reference 14 to similar CSF/vascular 
compliances as: 

Cab (P ab)=C A f(P ab) ^^=0. 02095 23 ml/mmHg 

C B APb v) = Cfv(Pb v) C^ r =0.688845 ml/mmHg. 

When calculating extra-cranial compliances it is advan- 
tageous to first determine the total compliance of each 
compartment. This quantity will be denoted by C z with the 
subscript indicating the compartment. This type of compli- 
ance can be described by the relation 

Total Compartmental Compliance = (78) 

Increase in Compartmental Volume 
Increase in Compartmental Pressure 


c gx= c x ml/mmHg (84) 

C^MaxCC^-^+C^^lml/mmHg (85) 

Here, the Max function is used to ensure that no negative 
40 compliances occur in the model. For the particular choices 
of total compliances here, this function is not needed, but it 
is included above as it may be required if total tissue 
compliances are calculated to be much smaller. 

Total blood volume (V Woorf ) is estimated at 5000 ml and 
45 5600 ml, where blood plasma makes up about 60% of this. 
The relative volumes of systemic veins, arteries and capil- 


laries are: 

V 5a =0.16 V b i ood systemic artery volume (86) 

50 V 5v =0.64 V blood systemic vein volume (87) 

V 5C =0.04 V b i ood systemic capillary volume (88) 

V^=0.09 V blood pulmonary system volume (89) 

Vhearr^-W ^ blood heart-blood volume (90) 


These values are based on estimates that the systemic 
arteries comprise 13% of the blood volume while the 
systemic arterioles and capillaries comprise 7%. In the 
60 present model, arteriole volume is lumped with the arteries, 
and hence 3% of the arteriole/capillary space is shifted here 
into the artery compartment resulting in the percentages 
above. It should also be noted that these percentages imply 
that systemic venous volume is four times that of systemic 
65 artery volume. 

As a next step, arterial and venous volumes must be 
allocated between the intracranial, lower, and central com- 
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partments. Estimates for this distribution are given below. It 
should be emphasized that these are percentages of the 
systemic blood and exclude the blood volumes in the heart 
and pulmonary system. 

pV cra „=0.07 percentage of systemic blood, intracra- 
nial region (91) 


pV /owr =0.40 percentage of systemic blood, lower 
region (92) 

pVcenfr-a/ = 0-53 percentage of systemic blood, central 
region 


10 


(93) 

The above percentages lead to the following volumes for 
the extra-cranial, vascular compartments: 


15 


v z= P viower V s a = l ower a^ery volume 

(94) 

V / =pV cew/ra/ V sa =central artery volume 

(95) 

Vj^pV /ower V sv =lower vein volume 

(96) 20 

V o=p V c<? „ Zr<3/ V J v =central vein volume 

(97) 

V Z) =pV /owr V s flower capillary volume 

(98) 

V 7 =pV ( , c „, ra/ V s<? =central capillary volume. 

25 

(99) 


The volumes of the tissue and rest of body compartments 
are now calculated based on interstitial fluid volume and 
intracellular fluid volume. Total interstitial fluid volume is 
estimated in the prior art to be 12 L, while total intracellular 30 
fluid volume is estimated to be 28 L. Considering that 
approximately 2 L of the intracellular fluid is contained in 
the blood cells, the following volume assignments can be 
made: 

35 

V^-12000 ml=interstitial fluid volume (100) 


V intra =26000 ml=intracellular fluid volume (101) 

With these volume assignments and the percentages in 
equations (92)-(93), estimates for lower tissue and rest-of- 
body volumes become: 

V G =P VloyveA VinteA E^Hcwer tissue volume 

V Y =PKentral(Vi n ter+Vi n tra)+Vpulrn=I * St of body VOl- 
ume 

Estimations for the pres sure- strain modulus are given by 
Milnor and Nichols and O’Rourke for various branches of 
the systemic arteries. From this data, the distensibility of the 
various vessels can be calculated. Of interest here are the 
calculated distensibilities for the ascending aorta and the 
femoral artery, represented in the model by D 7 and D z . 
Milnor’s data suggests that D 7 is 0.0036225 mmHg" 1 and D z 
is 0.00101 mmHg -1 while Nichols and O’Rourke data 
suggests D 7 is 0.00320141 mmHg -1 and D z is 0.00237912 55 
mmHg -1 . 

Averaging these values results in the assignments 

D/=0.00341196 (104) 60 

D^O.00169456. (105) 

Since it is known that the central veins are eight times as 
distensible as central arteries, D 0 is calculated accordingly 
by 65 


( 102 ) 
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(103) 


D 0 — 8D 7 


(106) 
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There are some suggestions in the prior art that the upper 
venous distensibility may be larger that of the lower leg, but 
these results are inconclusive. Therefore, lower venous 
distensibility is set to 

DjrDo. (107) 

Consistent with prior art findings, the capillary compart- 
ments C, J, and D are considered non-defoimable. 

The distensibility of the lower tissue compartment is 
assigned a weighted average of the interfacing vein and 
artery distensibilities. Since the veins have four times the 
volume of the arteries, these weights are used in 


D? I - 4 D y 

D g = — - = 0.0221754 (ml/mmHg)/ml 


(108) 


Finally, the distensibility of the rest-of-body compartment 
(Y) is also a weighted average based on its composition: 


P V «n, raider + V iam )D C + ( j D 0 + (5 V pu Jf>, 


(109) 


= 0.0221838 (ml/mmHg)/ml 


The distensibilities of the extra-cranial compartments 
have now been determined. Combined with the volumes of 
the previous section, each extra-cranial compartment now 
has an associated total compliance equal to the product of its 
associated volume and distensibility. 

A lymphatic system is also included in the present model. 
This system allows interstitial fluid in the tissue compart- 
ments of the central and lower regions of the body to flow 
directly into the central venous compartment. Lymphatic 
flows are thus able to balance filtration in these regions and 
establish an equilibrium state. Lymphatic flow Is not intro- 
duced in the intracranial region as the brain appears to lack 
lymphatic vessels. However, drainage of interstitial fluid via 
the Virchow-Robins Spaces is accommodated by pathways 
between the Brain and CSF compartments. 

The lymphatic system is represented in the present model 
through the flows Q ro and Q go . These flows provide path- 
ways through which interstitial fluid can move directly from 
the interstitium to the central venous system. Unlike the 
flows between adjacent compartments, which are driven by 
pressure differences, the lymphatic flow is governed almost 
exclusively by the interstitial fluid pressure. This flow 
increases 20 fold at interstitial fluid pressures near 4 mmHg 
from the corresponding flow at an interstitial pressure of -6 
mmHg. Below this interstitial fluid pressure value, lym- 
phatic flow ceases. The relative lymphatic curve is modeled 
by the logistic expression 


Plymph( x ) — 


M 

1 + (M - l)c-4*+6) 


where M = 20 and r = 1 


(110) 


This function, depicted in FIG. 4, defines the relative 
lymph flow in terms of interstitial pressure. The large dot in 
FIG. 4 denotes the mean pressure of interstitial fluid in the 
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model. Conversion of this relationship to the actual lym- 
phatic flows Q go and Q zo is accomplished by defining 


_ P lymph {Pg)QgO 

Qc0(t)= P^(Pc) ’ 

(HD 

^ ^ P lymph ( Py )Q V0 

firoW= P^CPr) ■ 

(112) 


In this manner, the mean flow is maintained at the mean 
pressure and can increase or decrease accordingly. Based on 
equation (110), Q GO can increase from a mean flow of about 
0.24 ml/min to 4.7 ml/min while Q yo can increase from 1 .76 15 
ml/min to 35 ml/min. 

Blood supply and CSF production is well regulated in the 
human intracranial system. This auto -regulation is achieved 
by the dilation and constriction of the blood vessels induced _ 
by factors such as carbon dioxide concentration, hydrogen 
ion concentration, oxygen concentration and intracranial 
pressures. The model incorporates such auto-regulation in 
the form of pressure sensitive fluidities depicting either 
vessel constriction in the form of a reduced fluidity or 25 
dilation in form of increased fluidity. 

In all, the present model contains 16 compliance param- 
eters, of which eight are variable and depend on pressure 
differences. There are also 23 resistance parameters. Four 
intracranial resistances are pressure-dependent and auto- 30 
regulate flow to both the cerebral capillary bed and the 
choroid plexus as well as the production of CSF in the 
choroid plexus. CSF production is known to be nearly 
constant over a wide range of pressures. However, the 
production 35 of CSF decreases when perfusion pressure 35 
decreases to low values, and a mechanism that ramps down 
CSF production when perfusion pressure falls below 55 
mmHg is also included in the model. 

Blood flow into the cerebral capillaries is auto-regulated 4Q 
by a pressure dependent fluidity Z AC defined by 


Z A c 


F(Pa ) 
Pa ~Pc 


(113) 
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where F(P^) is a logistic expression defined by 


P(Pa) = Q ac 


l+(M-l)c^- p ^ 


with M = 1.00001 and r — 0.25. 


(114) 


50 


This implies that Q^ C =Z^ C (P^ -P C )=F (P^) and cerebral 55 
blood flow is thus determined almost entirely by intracranial 
artery pressure and does not diminish significantly until this 
pressure drops below 50 mmHg. The relationship between 
Q^ c and P A is displayed in FIG. 5. 

Blood flow into the choroid plexus is regulated by a 60 
pressure-difference dependent fluidity Z AP defined by 




(115) 
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where the multiplier G(P per/ ) in equation (115) is a function 
of the perfusion pressure P^^P^-Pp that is unity when 
Pperf exceeds 55 mmHg and then falls linearly to zero below 
this value. This implies that 

Q A P=Z A AP A -Pp)=Q A pG(P per j). (116) 

The multiplier G has been included in equation (115) to 
model the fact that CSF production in the choroid plexus 
scales with blood flow, remains nearly constant over a wide 
range of pressures, and is decreased as a function of the 
magnitude of the perfusion pressure when P perf falls 20 
below 50 to 60 mmHg. From equation (116), blood flow into 
the choroid plexus remains constant unless perfusion pres- 
sure falls to a low level. 

Over the wide range of normal pressures, the production 
of CSF is auto -regulated at the venous side of the choroid 
plexus capillaries by the pressure dependent fluidity Z PV 
defined by 


Zpv = 


Q AP -G{P petf )-Z PF -{P P -P F ) 

Pp-Pv 


(117) 


When V per j^55, this expression for Z PV will maintain a 
constant pressure difference between the choroid plexus and 
the ventricular CSF. Substituting Z PV with G=1 into the 
governing equation for the choroid plexus, equation (15) 
reduces to 


d Ppp 
dt 


(118) 


Since the compliance C PF must be non-zero to account for 
the known ability of the Choroid Plexus to transmit pressure 
pulsations to the ventricular CSF, the governing equation for 
the choroid plexus compartment becomes simply 


dPpF =() (119) 

dt 


This implies a constant pressure difference between the 
choroid plexus and ventricular CSF is maintained by equa- 
tion (117) for V perf ^55 mmHg. Therefore 

Q PF =Zpp{P P -P p) =Z PF (P P -P p)=Q P p. (120) 

Since for pressures in the normal range, CSF production 
in the choroid plexus is proportional to P PF , constant CSF 
production from the choroid plexus is thus achieved. Equa- 
tion (119) also eliminates the need to estimate C PF in this 
model as occurrences of this parameter are always multi- 
plied by either 6P FF Jd\ or dPpp/dt. 

Using the above results, it can now be demonstrated how 
Z pv autoregulates CSF production for Pp ery ^55 mmHg. 
Substituting Q PF for Z pf (P p -P f ) in equation (117) and 
noting that Q^p-Qpp=Qp F results in the equalities: 


Zpv 


Qpy 

Pp-Pv 


Qpy 


(121) 


65 


(Pp - P F ) + (P F - Py) 
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-continued 
Qpv 

Q PF / z PF + ( Pp -Py) 


The last term in this expression reveals the relationship 
between Z pv and P^. Physiologically, Z PV should decrease 
with increasing CSF pressure (P^) causing an increase in 
choroid plexus pressure and maintaining constant CSF pro- 
duction across Z PF- It is quite clear from the last equality in 
equation (121) that this is indeed the case. These autoregu- 
latory mechanisms may be substituted into (13) to (16) and 
(18) to obtain the governing equations for compartments A, 
C, P, V, and F. 

The present model also includes a group of regulatory 
mechanisms associated with the sympathetic nervous sys- 
tem (SNS). Two variable resistances in the central and lower 
regions provide for SNS regulation of arterial pressure 
through constriction of the arterioles. Included in these 
variable resistances is the dramatic “last ditch stand” trig- 
gered by the SNS when arterial pressure falls to a level 
where cerebral blood flow is affected. The far less extreme 
SNS regulation of arterial pressure through a constriction of 
the large vascular vessels in the central body is also repre- 
sented in the model. Active interfaces are placed between the 
central rest-of-body compartment and the central artery and 
vein compartments. When arterial pressure falls, forcing 
terms in the governing equations for compartments Y, I, and 
O force the volume cups at the active Y-I and Y-O interfaces 
into the vascular compartments, providing the regulatory 
constriction. An additional SNS mechanism in the model 
regulates central arterial pressure by increasing the number 
of heartbeats per minute if arterial pressure falls. 

The sympathetic nervous system (SNS) is associated with 
reflex mechanisms that act to maintain levels of arterial 
pressure and cardiac output if arterial pressure falls. This 
section deals with the portions of the model that represent 
the capacity of the SNS to rapidly increase arterial pressure 
by constricting both the arterioles and the large vascular 
vessels. The model’s regulation of cardiac output by the 
SNS will be developed below. 

A first SNS pressure mechanism involves constriction of 
the large vessels of the circulation in the central body. If 
arterial pressure falls, the large vessels of the circulation in 
the central body, especially the veins, strongly constrict to 
cause a rapid increase in arterial pressure. This mechanism 
is included in the current model by placing active interfaces 
between the central compartment Y and the central vascular 
compartments I and O. Forcing terms in the conservation 
equation for compartment Y now force the volume cups at 
the interfaces to push into the I and O compartments when 
arterial pressure diminishes, modeling the SNS large vessel 
constriction mechanism. 

The conservation equation in compartment Y is described 
by 


dVy dV YO dVyi dVy M (122) 

-dT = — + — + — =Qjr - Qm 


where the second equality determines the equation. Here, the 
volume change dV^dt between compartment Y and the 
external environment M, as in equation (3), is simply 
proportional to the change in pressure difference dP^^/dt. 
However, dV FO /dt and dV^dt involve both changes in the 
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pressure differences dP FO /dt and dP F7 /dt, respectively, and a 
forcing term describing the active compliance between Y 
and O and between Y and I, respectively. In particular, 


5 


d Vyo d Pyo 

~dT~ ~ Fyo + Ho(0 with F y 0 {t) = -24 A 
and 


(123) 


10 


d Vyi d Pyi 

~dT = Cy, ~dT + Fyi{t) with Fyi(t) = ~ 8p/ ' 


(124) 


The governing differential equation in compartment Y is 
now obtained by replacing dV FO /dt and dV 17 /dt in equation 
15 (7) by the expressions (123) and (124). Similarly, for the 
governing equations in compartments O and I, dV FO /dt is 
replaced in equation (8) by expression (123) and dV^/dt is 
replaced in equation (5) by expression (124), respectively. 

A second SNS pressure regulation mechanism involves 
20 constriction of the arterioles in most parts of the body 
(excluding the brain and heart) when central artery pressure 
drops. This causes an increase in arterial pressure through an 
increase in the total peripheral resistance. When the arterial 
pressure drop is severe enough to compromise blood flow to 
25 the brain, this regulatory response is extreme. 

To model this portion of the SNS regulatory response, two 
multipliers are defined for the artery-capillary fluidities Z u 
and Z IZD . These variable resistances are of the form 

30 

Qn (125) 

Z u = =^4- ■ SNSzi (Pi) ■ SNSz(Qac) 

Pi-Pj 

and 

35 Qzd (126) 

Zzd = _ ^ ■ SNSzi (P/) ■ SNSz(Q ac ) 

Pz - Pd 


The first multiplier, 
40 


SNSzi(Pi)=j, 


(127) 


45 

is a function of central body artery pressure that increases 
resistance if P 7 falls below its mean level. The second 
multiplier, which is a function of the cerebral blood flow 
Q 4C , is defined by 

50 


55 


SNSz(Q ac ) = 


M 

1 +(M-l)c^ A c-e A c) 


where 


M = l.l and r = 0.02. 


(128) 


This multiplier models the last-ditch stand when cerebral 
blood flow is significantly reduced. SNSz remains close to 
unity (not activated) until arterial pressure drops to a level 
60 where Q AC is affected. At this point, SNSz drops sharply, 
dramatically increasing the resistances R u and R zz> 

The major source of forcing in the present whole-body 
model comes from the heart and involves the cardiac output 
and cardiac uptake Q OH . Because it is known that all 
65 extra blood that flows into the heart is automatically pumped 
without delay into the aorta and flows again through the 
circulation, the cardiac output will be set equal to cardiac 
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uptake, so that Q^Q OH It is farther known that the normal 
heart, functioning without any special stimulation, can pump 
an amount of venous return up to 2.5 times the normal 
venous return before the heart becomes the limiting factor. 
Therefore, a venous return function R is defined that incor- 5 
porates all of the flow into the central venous compartment 

R=Qyo+Qgo+Zjo(P^P o)+Z to (P t~P o)+Z so (P s~ 

Po)+Z X o(PjtPO) (129) 

The cardiac uptake Q OH is now defined in terms of this 
venous return and two regulatory multipliers SNSo and OVP 
by 

Q OH =SNS o (P',)-0rP(P o )-R (130) 

where P ' 7 denotes the time derivative of the arterial pressure 15 
function P 7 . Since Q hi = Q 0 h’ cardiac output is now based on 
the venous return through equation (130). 

The multiplier SNSo in equation (130) is associated with 
cardiac regulation by the sympathetic nervous system 
(SNS). It models an SNS mechanism that increases or 20 
decreases the heart rate in response to a change in arterial 
pressure. A decrease in arterial pressure can markedly 
increase heart activity, with the heart rate increasing to as 
high as 160-180 bpm from the normal value of 72 bpm. A 
linear model for the SNSo multiplier as a function of the 25 
instantaneous time derivative P' 7 is 


w P '1 ( 131 ) 

SNSoiP',) = 1 - jL. 

This relationship produces an increase in heart rate to 
about twice the normal level as the result of a pressure drop 
of 20 mmHg. Note that SNSo (0)=1, so that this regulatory 
mechanism is only active when P 7 is changing. 

The OVP function in equation (130) insures that if venous 
pressure drops, then so does cardiac uptake. It is defined by 


OVP(Po) = 


M 

1 + (M - 1 )e r<p o- p o) 


where M = 2.5 and r — 0.5. 


(132) 40 


Since OVP(P 0 )=1, this regulatory mechanism is not active 
when central venous pressure remains at its mean value. 

Using the above cardiac forcing terms leads to a system 
of equations that describe mean pressures in the sense that 
the oscillations that occur about a mean pressure value are 50 
removed and the mean pressure may be considered a tem- 
poral average over one cardiac cycle. However, when trying 
to resolve circulatory pressure pulsations caused by cardiac 
output, a more instantaneous description of cardiac output is 
necessary. Again in this case, cardiac output is set equal to 55 
cardiac uptake. 

Referring now to FIG. 6, it is contemplated that math- 
ematical model 10 of the present invention will operate in a 
computing environment including a computer system 20. 
Computer system 20 includes a computer central processing 60 
unit (CPU) 22, a computer memory 24, and input/output 
devices 26. Mathematical model 10 is typically imple- 
mented with or embodied in a computer program 28 which, 
when executed by computing resources within computer 
system 20, provide the functionality of the present invention. 65 
Typically, computer program 28 resides in computer 
memory 24 of an individual client computer system. Of 


course, computer program 28 may reside in the memory of 
a local or wide area network server or in the memory of an 
equipment computer processor. In one embodiment, com- 
puter programs 28 include a body compartment module 29, 
a flow equation module 30, and an equation solver module 
31. Body compartment module 29 includes instructions for 
dividing the body into a plurality of compartments and a 
heart pump, each compartment representing a portion of the 
body. Flow equation module 30 includes instructions for 
deriving a plurality of differential flow equations, each of 
which correspond to one of the compartments. Equation 
solver module 31 includes instructions for solving the plu- 
rality of differential flow equations. The instructions within 
body compartment module 29, flow equation module 30, 
and equation solver module 31 are executed within com- 
puter programs 28 to simulate the pressure dynamics of the 
intracranial system. Input/output devices 26 typically 
include a storage device 36, such as a hard disk drive, a 
keyboard 38, a pointing device 40, i.e., a mouse, and a 
display device 42, such as a monitor. It is contemplated that 
data may be manually input to mathematical model 10 via 
input/output devices 26 or automatically input via patient 
monitoring equipment or the like. 

To validate the present whole-body model for intracranial 
pressure dynamics, several types of simulations assuming 
normal physiology were carried out and the results com- 
pared to clinical observations. In the first type of simulation, 
a pulsatile form for cardiac output was introduced as forcing 
so that the model’s predicted pressure responses could be 
compared to clinically measured pressure pulsations in 
healthy humans. The pulsatile cardiac output function with 
parameter values n=ll, <|)=0.27882, and cr=5.958 was used 
as the cardiac forcing function (Q HI ) in the model’s gov- 
erning equations. This function is depicted in FIG. 7 A and 
represents the cardiac output over two cycles. A mean value 
for central artery pressure of 96 mmHg was prescribed as an 
initial condition for the first simulation. In the second type 
of simulation, a constant flow infusion test was used to 
validate the model’s representations for variable intracranial 
compliances and intracranial fluid dynamics. In these simu- 
lations, the governing equation for the extra-ventricular CSF 
compartment (T) was augmented by adding a constant 
infusion term to model the clinical infusion of mock CSF 
into the lower lumbar space. 

With the lower body region explicitly represented by 
separate compartments in the present model, the effect of 
orthostatic forces on cerebral blood flow as the body 
changes position can be considered. Two additional types of 
simulations were now run to validate the modeling of the 
regulatory mechanisms associated with the systemic ner- 
vous system. The normal physiology value for the resistance 
R xo was increased twofold and R /z was decreased by one 
half to simulate a change in body position from lying down 
to standing up, and the behavior of the cerebral blood flow 
Q^ c was then examined. These resistance changes were 
made instantaneously and the modeled effect on cerebral 
blood flow was determined for the next 30 seconds. In the 
first simulation of this type, the model equations included all 
of the SNS reflexes described above. In the second simula- 
tion, the SNS terms were removed from the governing 
equations. 

To examine the potential predictive capabilities of the 
model in pathological conditions, simulations were run 
where all cardiac forcing in the model was suddenly stopped 
to simulate cardiac arrest, i.e., and Q OH were suddenly 

set equal to zero. A blood volume of 5600 ml was assumed 
in these simulations, and the response of pressures in the 
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model’s circulatory compartments was determined. This 
behavior, and the predicted final circulatory compartmental 
pressure values, were then compared to clinical results 
associated with cardiac arrest and the filling pressure of the 
circulation. 

As a second example of pathology, simulations of hem- 
orrhagic shock were carried out. Hemorrhage was modeled 
by the inclusion of an outflow path, denoted Q XM , from 
lower venous compartment X into ambient atmosphere 
compartment M. This flow was calculated so as to achieve 
a 45% loss in blood volume at the end of the simulation. The 
percent changes in central artery pressure, cardiac output, 
and cerebral blood flow were then calculated with respect to 
percent blood loss. 

In all of these simulations, the model’s system of differ- 
ential equations was solved numerically using a typical 
symbolic mathematical software package employing maxi- 
mum accuracy settings. 

The usual systolic and diastolic values of pulsatile central 
artery pressure are perhaps the best-known values in human 
physiology. FIG. 7B shows the behavior of central artery 
pressure predicted by the present model in response to the 
pulsatile cardiac forcing developed in the prior art. The 
model’s predicted response to a mean value of 96 mmHg, 
prescribed as an initial condition, reproduces a peak systolic 
pressure of 120 mmHg and a diastolic pressure of 80 mmHg. 
FIG. 7B also depicts the intracranial artery pressure 
response. Here, the predicted systolic and diastolic pressures 
are “100 over 65.” These values are typical of those in the 
small arteries. 

Because even the most complex mathematical model 
must be based on assumptions and simplifications of actual 
physiology, model validation is an essential step in the 
development process. In the present case, after calibration of 
parameters associated with healthy human physiology, the 
model was used in two types of simulations and the results 
compared to physical data. In the first type of simulation, the 
response of compartmental pressures to the realistic pulsatile 
cardiac output given in FIG. 7 A was determined. As shown 
in FIG. 7B the predicted response of central arterial pressure 
is the typical “120 over 80” blood pressure readings 
expected in healthy humans. It is worth noting that a mild 
incisura (dicrotic notch) is discemable in this pressure curve. 
This notch is typical in central artery pressure readings, 
although is it usually more pronounced than in FIG. 7B. In 
particular, the present model cannot capture reflected waves 
which have a tendency to cause a second (usually higher) 
systolic peak and hence a more prominent inflection point. 

The response of the intracranial arteries is also shown in 
FIG. 7B. This response agrees well with prior art estimates 
for pressure pulses in the small arteries. The pressure 
responses of other intracranial compartments were also 
within expected ranges. 

Constant flow infusion tests were also simulated using the 
present model. In these experiments, mock CSF was infused 
at a constant rate into the lower lumbar space. The pressure 
of this space was then measured and associated with a 
calculated total CSF volume change. This resulted in deter- 
mination of a curve known as the global pressure-volume 
relation. The slope of this curve describes the elastance of 
the entire CSF space, including extracranial portions. The 
inverse of the elastance is the more well-know compliance. 

The typical clinical global CSF pres sure- volume relation, 
except at extreme pressures, is an S-shaped curve of logistic 
type. It has a lower pressure plateau near resting pressure, 
defined as the pressure where CSF production in the system 
is just balanced by CSF absorption by the venous system. 
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This region of small slope (large compliance) is due to the 
ability of the system to easily accommodate increases in the 
volume of the CSF space at these relatively low pressures 
through the compression of the venous system (excluding 
5 the venus-sinus veins). As additional CSF volume is added 
to the system and pressures increase, this capacity for 
adjustment diminishes as there will be less venous blood 
available to eject to further compress the veins. Thus, with 
increasing pressures, the pres sure- volume curve steepens 
to indicating a reduction in the compliance of the system. The 
slope of the pres sure- volume relationship continues to 
increase for larger infusion volumes until the resulting CSF 
pressures are high enough that the intracranial arteries can 
begin to be compressed. At this point, some additional 
15 compliance enters the system. There is a point of inflection 
in the curve followed by a region of decreasing slope that 
leads to an upper pressure plateau at the diastolic pressure of 
the intracranial arteries. Once additional volume increases 
beyond this point increase CSF pressures to the systolic 
20 pressure of the intracranial arteries, there are no additional 
mechanisms available to buffer additional volume increases, 
and the compliance of the CSF system falls to zero. 

The pressure-volume relation obtained by using the 
present model to simulate CSF infusion tests is given in FIG. 
25 8 A. The pressure in this figure is the pressure P r of the 
extraventricular CSF compartment since this compartment 
contains the lower lumber portion of CSF space where 
clinical pressure recordings were made. The volume change 
in FIG. 8 A represents the net change of both ventricular and 
30 extra-ventricular CSF volumes combined, as this is the 
quantity calculated in the experiments. The logistic -like 
shape of the predicted pressure- volume relationship is main- 
tained until CSF pressures reach high values comparable to 
the systolic blood pressure. At this point, the curve turns 
35 upward and compliance rapidly decreases. 

The same infusion simulations that led to the pressure- 
volume relationship in FIG. 8 A may also be used to calculate 
a relationship between CSF pressure increases and CSF 
absorption. The slope of this relationship is known as the 
40 conductance of CSF outflow and is denoted by C out . In 
model simulations, this value is easily calculated by 

ZTs(PT(t)-Ps(t))+ ( 133 ) 

45 _ n _ z TO (p T ( t ) - pm) - (q ts + e ro ) 


50 where the superscript m denotes the model calculation and 
t denotes the time in minutes from the start of the simulated 
infusion. FIG. 5B illustrates the predicted conductance of 
CSF outflow throughout ten minutes of an infusion simula- 
tion. 

55 When C out has been has been studied experimentally, a 
linear relationship between CSF pressure increase and CSF 
absorption is observed. In prior art experiments, a mean 
value for C out of 0.11 (mlmin)/mmHg is given for a sample 
of eight healthy volunteers, and it is stated that a value 
60 greater than 0.10 is probably normal. The values of the 
conductance of CSF outflow calculated from the present 
simulations are shown in FIG. 8 B. These values change with 
time, but stay within 0.004 units of the mean value of 0.11 
(ml/min)/mmHg observed in the prior art. They are also 
65 greater than 0.10 as previously suggested. 

Furthermore, the calculated temporal variation of C out m is 
sufficiently small that that the relationship between CSF 
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pressure increase and CSF absorption might easily be cat- 
egorized as linear on the basis of clinical data. 

In the current model calibrations, the pressures, flows, and 
hence resistances are determined from data associated with 
the body being in the horizontal (supine) position. Gravita- 
tional influences on the circulation will change when the 
body changes position. Upon standing up, the blood flow 
into the lower body is aided by gravity, but the venous return 
from the lower body is hindered by it. As the lower body in 
the present model is represented by separate compartments, 
the effect of a positional change can be considered. To 
simulate the gravity-induced changes associated with stand- 
ing up, the resistance into the lower arteries (R /z ) may be 
decreased by one half while the resistance from the lower 
veins (Rj^o) * s doubled. This results in pressure and volume 
increases of the lower arteries and veins. In the current 
simulations, these resistance changes were made instanta- 
neously and the modeled effect on cerebral blood flow was 
determined for the next 30 seconds. In the first simulation, 
all sympathetic nervous system (SNS) reflexes are intact and 
in the second 39 these are removed. As can be seen in FIG. 
9, with all SNS reflexes activated, cerebral blood flow is 
maintained at 100% of its original value. 

Referring now to FIG. 9, the results of modifying the 
resistances between the central and lower body in the model 
to simulate a change in body position from lying down to 
standing up are illustrated. The solid curve in FIG. 9 
indicates the response of cerebral blood flow to this change 
in position with the SNS reflexes included in the model 
equations. The simulation predicts that with all SNS reflexes 
activated, cerebral blood flow is maintained at 100% by the 
model through the change in body position. The dashed 
curve in FIG. 9 shows the predicted behavior of cerebral 
blood flow in the absence of regulation by the SNS. 

When the SNS reflexes are removed from the model 
equations, cerebral blood flow drops to about 94% of its 
original value due to the positional change. This predicted 
decrease indicates the important role that the SNS regulatory 
mechanisms play in the maintenance of cerebral blood flow 
in the model. The decrease also confirms that loss of the SNS 
reflexes can be a factor in orthostatic intolerance. Syncope, 
or fainting, is a common event with many different causes. 
Syncope can be due to sudden vasodilatation (vasodepressor 
or “vasovagal” syncope) as well as postural hypotension, in 
which the normal vasoconstrictive reflex response to a 
transiently decreased cardiac output on standing is not 
sufficiently active. Postural hypotension can occur as the 
result of drugs, venous disease, sympathectomy, hypov- 
olemia, peripheral neuropathy, in addition to degeneration of 
the sympathetic nervous system (primary autonomic insuf- 
ficiency, or idiopathic orthostatic hypotension). Common to 
all causes of syncope is a decrease in cerebral blood flow to 
under 30 ml per 1 00 g brain tissue per minute from the usual 
50-55 ml, about 55% of normal. This is well above the 
threshold for loss of electrical function, at 30% of normal, 
but is enough to cause a transient loss of consciousness. The 
calculated percentage decrease in cerebral blood flow in the 
current simulation with the SNS terms in the model equa- 
tions omitted does not approach the levels associated with 
fainting due to a rapid change in position as the result of 
standing up too quickly, even with inactive SNS reflexes. 
However, cerebrovascular autoregulation in the current 
simulation remains uncompromised and will act to maintain 
cerebral blood flow despite sudden decreases in arterial 
pressure. 

To demonstrate the potential predictive capabilities of the 
model in pathology, two situations were simulated. In the 
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first of these, a simplified representation of cardiac arrest 
was created by suddenly terminating all cardiac forcing, i.e., 
setting Qtf/=Qoiy = 0, at a specified time in the course of the 
simulation. In prior research by Gutyon, it is noted that 
5 “When heart pumping is stopped . . . the flow of blood 
everywhere in the circulation ceases a few seconds later. 
Without blood flow, the pressures everywhere in the circu- 
lation become equal after a minute or so. This equilibrated 
pressure level is called the mean circulatory filling 
to pressure ... at a volume 40 of 5000 ml, the filling pressure 
is the normal value of 7 mmHg.” Predicted results from the 
simulation were consistent with this statement. When all 
cardiac forcing was stopped in the simulation, all circulatory 
pressures tend to 7.5 mmHg after about 45 seconds. The 
15 small difference in the predicted and cited filling pressures 
may be partially due to the fact that total blood volume in the 
simulation was taken to be 5600 ml rather than the 5000 ml. 

Referring now to FIG. 10, the second pathological situ- 
ation considered as a test of the model’s capabilities 
20 involves the effects of shock caused by hypovolemia. In 
these simulations, a flow term was introduced into the 
equation for compartment X to model a hemorrhage from 
the lower body. Prior art research by Guyton gives a 
clinically-derived graph that, as in FIG. 10, depicts the 
25 percentage drop in central artery pressure (P 7 ) and cardiac 
output (Qjjj) with respect to percent blood loss. Consistent 
with these clinical results, the modeled response shows very 
stable pressure and cardiac output over the first 10% blood 
loss, and the relative arterial pressure stays above relative 
30 output during the entire course of the hemorrhage. At 30% 
blood loss, the relative drops in pressure and cardiac output 
are also in close agreement with the clinical results. During 
the course of the simulated hemorrhage, the central venous 
compartment contracts to less than 50% of its original 
35 volume, indicating that the active interface (SNSc) in the 
model is causing the central veins to strongly constrict in 
order to maintain blood flow back to the heart. Beyond 30% 
blood loss, a second arterial plateau is noticed in FIG. 10. 
This important feature, which is also noted in the clinical 
40 results, is due to the “last-ditch stand” reflex to maintain 
blood supply to the brain provided by the regulatory mul- 
tiplier SNSz (Q^ c ) in equation (130). 

Referring now to FIG. 11, the predicted percentage drop 
in cerebral blood flow (Q^ c ) with respect to percent blood 
45 loss during the simulated hemorrhage is illustrated. The 
solid curve represents results when both SNSz mechanisms 
are intact, and the dashed curve represents results when both 
SNSz mechanisms are disabled. Both simulations predict 
that a 25% blood loss can be tolerated without a significant 
50 reduction in cerebral blood flow. Above this degree of blood 
loss, with the SNS mechanisms intact, the blood supply to 
the brain remains above 75% of its original value. However, 
without these SNS regulatory mechanisms included in the 
model’s equations, cerebral blood flow drops quickly to 
55 20% its original value. 

The role of the SNS regulatory mechanisms in maintain- 
ing adequate central artery pressure and cardiac output has 
been noted above. However, the need to include the arteriole 
constrictive reflex SNSz in a model of intracranial pressure 
60 dynamics when simulating pathology is also dramatically 
illustrated by the behavior of cerebral blood flow as a 
function of blood loss in the hemorrhage simulations. Due to 
the intracranial autoregulatory mechanism Z AC in the model, 
a 25% blood loss can be tolerated without significant reduc- 
65 tion in cerebral blood flow. Beyond this point, arteriole 
constriction is essential to further maintain adequate levels 
of cerebral blood flow. As indicated in FIG. 11, predicted 
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cerebral blood flow Q^ c remains above 80 % of its original 
value from a 25 % blood loss through a 40 % blood loss. By 
contrast, in an analogous simulation where the arteriole 
constrictive response was disabled, Q AC falls to approxi- 
mately 55 % of its original value when a 40 % blood loss has 5 
occurred. Between a 40 % and 45 % blood loss, Q^ c remains 
above 75% of its original value with the arteriole constric- 
tive response intact, but falls rapidly from 55 % to approxi- 
mately 20% of its original value with the arteriole constric- 
tive response disabled. 10 

Most attempts to study intracranial pressure using 
lumped-parameter models have adopted the classical 
“Kellie- Monro Doctrine,” which considers the intracranial 
space to be a closed system that is confined within the 
nearly-rigid skull, conserves mass, and has equal inflow and 15 
outflow. The present work revokes this Doctrine and devel- 
ops a mathematical model for the dynamics of intracranial 
pressures, volumes, and flows that embeds the intracranial 
system in extensive whole-body physiology. The new model 
consistently introduces compartments representing the tis- 20 
sues and vasculature of the extradural portions of the body, 
including both the thoracic region and the lower extremities. 

In addition to vascular connections, a spinal subarachnoid 
cerebrospinal fluid (CSF) compartment bridges intracranial 
and extracranial physiology allowing explicit buffering of 25 
intracranial pressure fluctuations by the spinal theca. The 
model contains cerebrovascular autoregulation, regulation 
of systemic vascular pressures by the sympathetic nervous 
system, regulation of CSF production in the choroid plexus, 
a lymphatic system, colloid osmotic pressure effects, and 30 
realistic descriptions of cardiac output. 

Validation tests show that the model’s predictions are in 
agreement with experimental data for normal physiology. 
Additional simulations indicate that the present whole-body 
model appears to have significant potential predictive capa- 35 
bilities in situations involving pathological conditions. 

One embodiment of the present invention has been dis- 
closed and discussed herein, although it should be under- 
stood that the present invention is not limited to this (or any 
other) particular embodiment. On the contrary, the present 40 
invention is intended to cover all alternatives, modifications 
and equivalents that may be included within the spirit and 
scope of the appended claims. 

What is claimed is: 45 

1. A method of mathematically modeling pressure dynam- 
ics of a body’s intracranial system comprising the steps of: 
dividing the body into a plurality of compartments and a 
representation of a heart pump, said plurality of com- 5Q 
partments including at least one intracranial compart- 
ment and at least one extracranial compartment, each of 
said plurality of compartments representing a portion of 
the body, said representation of a heart pump interact- 
ing with at least one of said plurality of compartments; 55 
deriving a plurality of differential equations, each of said 
plurality of differential equations governing a pressure 
dynamic of one of said plurality of compartments; 
defining an atmosphere compartment, said atmosphere 
compartment representing a space located outside the 60 
body and wherein at least one of said plurality of 
differential equations accounts for a pressure of said 
atmosphere compartment; 

solving said plurality of differential equations; and 
producing a modeled circulatory system dynamic of the 65 
body based on said solving of said plurality of 
differential equations. 
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2. A method according to claim 1, wherein at least one of 
said plurality of differential equations accounts for cere- 
brovascular autoregulation byte body’s sympathetic nervous 
system. 

3. A method according to claim 1, wherein said plurality 
of compartments include a plurality of vascular compart- 
ments. 

4. A method according to claim 3, wherein said plurality 
of vascular compartments include a plurality of intracranial 
compartments. 

5. A method according to claim 4, wherein said plurality 
of intracranial compartments represent at least one of the 
intracranial arteries, intracranial capillaries, choroids plexus 
capillaries, venous sinus jugular veins, and intracranial 
veins. 

6. A method according to claim 3, wherein said plurality 
of vascular compartments include a plurality of central body 
compartments. 

7. A method according to claim 6, wherein said plurality 
of central body compartments represent at least one of the 
central arteries, central capillaries, central veins, and extra- 
ventricular CSF. 

8. A method according to claim 3, wherein said plurality 
of vascular compartments include a plurality of lower body 
compartments. 

9. A method according to claim 8, wherein said plurality 
of lower body compartments represent at least one of the 
lower arteries, lower capillaries, and lower veins. 

10. A method according to claim 1, wherein said plurality 
of compartments include a plurality of non-vascular com- 
partments. 

11. A method according to claim 10, wherein said plural- 
ity of non-vascular compartments represent at least one of 
the lower tissue, brain, ventricular CSF, and extra-ventricu- 
lar CSF. 

12. A method according to claim 1, wherein said plurality 
of compartments include a rest of body compartment. 

13. A method according to claim 1, wherein at least one 
of said plurality of differential equations include a term 
representing pressure flow into and/or out of a correspond- 
ing one of said plurality of compartments. 

14. A method according to claim 1, wherein said plurality 
of differential equations include an equation simulating fluid 
filtration from capillaries into interstitial space. 

15. A method according to claim 1, wherein said plurality 
of differential equations include an expression simulating 
deformation of the membrane between adjacent compart- 
ments. 

16. A system for modeling an intracranial system com- 
prising: 

a body compartment module adapted to divide the body 
into a plurality of compartments and a representation of 
a heart pump, said plurality of compartments including 
at least one intracranial compartment and at least one 
extracranial compartment, each of said plurality of 
compartments representing a portion of the body; 

an atmosphere compartment, said atmosphere compart- 
ment representing a space located outside the body; 

an equation module adapted to derive a plurality of 
differential equations, each of said plurality of differ- 
ential equations governing a pressure dynamic of one 
of said plurality of compartments; 

an equation solver module adapted to solve said plurality 
of differential equations wherein at least one of said 
plurality of differential equations accounts for a pres- 
sure of said atmosphere compartment; and 
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an output device adapted to output a result from said 
equation solver as a model of a circulatory system 
dynamic of the body. 

17. A system according to claim 16, wherein at least one 
of said plurality of differential equations accounts for cere- 
brovascular autoregulation by the body’s sympathetic ner- 
vous system. 

18. A system according to claim 16, wherein said plurality 
of compartments include a plurality of vascular compart- 
ments. 

19. A system according to claim 18, wherein said plurality 
of vascular compartments include a plurality of intracranial 
compartments. 

20. A system according to claim 19, wherein said plurality 
of intracranial compartments represent at least one of the 
intracranial arteries, intracranial capillaries, choroids plexus 
capillaries, venous sinus jugular veins, and intracranial 
veins. 

21. A system according to claim 18, wherein said plurality 
of vascular compartments include a plurality of central body 
compartments. 

22. A system according to claim 21, wherein said plurality 
of central body compartments represent at least one of the 
central arteries, central capillaries, central veins, and extra- 
ventricular CSF. 

23. A system according to claim 18, wherein said plurality 
of vascular compartments include a plurality of lower body 
compartments. 

24. A system according to claim 23, wherein said plurality 
of lower body compartments represent at least one of the 
lower arteries, lower capillaries, and lower veins. 

25. A system according to claim 16, wherein said plurality 
of compartments include a plurality of non-vascular com- 
partments. 

26. A system according to claim 25, wherein said plurality 
of non-vascular compartments represent at least one of the 
lower tissue, brain, ventricular CSF, and extra-ventricular 
CSF. 

27. A system according to claim 16, wherein said plurality 
of compartments include a rest of body compartment. 

28. A system according to claim 16, wherein at least one 
of said plurality of differential equations include a term 
representing a pressure flow into and/or out of a correspond- 
ing one of said plurality of compartments. 

29. A system according to claim 16, wherein said plurality 
of differential equations include an equation simulating fluid 
filtration from capillaries into interstitial space. 

30. A system according to claim 16, wherein said plurality 
of differential equations include an expression simulating 
deformation of the membrane between adjacent compart- 
ments. 

31. A method of modeling pressure dynamics of an 
intracranial system comprising the steps of: 

dividing a body into a plurality of compartments and a 
representation of a heart pump, each of said plurality of 
compartments representing a portion of the body, said 
plurality of compartments including at least one intrac- 
ranial compartment and at least one extracranial com- 
partment, said representation of a heart pump interact- 
ing with at least one of said plurality of compartments; 

establishing an atmosphere compartment, said atmo- 
sphere compartment representing a space located out- 
side the body; 

deriving a plurality of differential equations, each of said 
plurality of differential equations governing a pressure 
dynamic of one of said plurality of compartments, 
wherein at least one of said differential equations 
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accounts for cerebrovascular autoregulation by a sym- 
pathetic nervous system, wherein at least one of said 
plurality of differential equations accounts for a pres- 
sure of said atmosphere compartment; 

5 solving said plurality of differential equations; and 

producing a modeled circulatory system dynamic of the 
body based on said solving of said plurality of differ- 
ential equations. 

32. A method of modeling pressure dynamics of an 
to intracranial system comprising the steps of: 

dividing a body into a plurality of compartments, each of 
said plurality of compartments representing a portion of 
the body, said plurality of compartments including at 
least one intracranial compartment and at least one 
15 extracranial compartment; 

establishing an atmosphere compartment, said atmo- 
sphere compartment representing a space located out- 
side the body; 

providing a means for representing a heart pump that 
20 interacts with at least one of said plurality of compart- 
ments; 

deriving a plurality of differential equations, each of said 
plurality of differential equations governing a pressure 
dynamic of one of said plurality of compartments, 
25 wherein at least one of said differential equations 
includes a means to account for cerebrovascular auto- 
regulation by a sympathetic nervous system, wherein at 
least one of said plurality of differential equations 
accounts for a pressure of said atmosphere compart - 
30 ment; and 

solving said plurality of differential equations; and 
producing a modeled circulatory system dynamic of the 
body based on said solving of said plurality of differ- 
ential equations. 

35 33. A method of modeling pressure dynamics of an 

intracranial system comprising the steps of: 
providing a means for dividing a body into a plurality of 
compartments and a representation of a heart pump, 
each of said plurality of compartments representing a 
40 portion of the body, said plurality of compartments 
including at least one intracranial compartment and at 
least one extracranial compartment, said representation 
of a heart pump interacting with at least one of said 
plurality of compartments; 

45 providing a means for establishing an atmosphere com- 
partment, said atmosphere compartment representing a 
space located outside the body; 
providing a means for deriving a plurality of differential 
equations, each of said plurality of differential equa- 
50 tions governing a pressure dynamic of one of said 
plurality of compartments and said atmosphere com- 
partment; and 

providing a means for solving said plurality of differential 
equations and producing a modeled circulatory system 
55 dynamic of the body based on said solving of said 
plurality of differential equations. 

34. A method of modeling pressure dynamics of a body’s 
intracranial system comprising the steps of: 

dividing the body into a plurality of compartments and a 
60 representation of a heart pump, each of said plurality of 
compartments representing a portion of the body, said 
plurality of compartments including at least one intrac- 
ranial compartment and at least one extracranial com- 
partment, said representation of a heart pump interact- 
65 ing with at least one of said plurality of compartments, 
wherein a plurality of said plurality of compartments 
are vascular and a plurality of said plurality of com- 
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partments are non-vascular, said vascular compart- 
ments including at least one of the intracranial arteries, 
intracranial capillaries, choroids plexus capillaries, 
venous sinus jugular veins, intracranial veins, central 
arteries, central capillaries, central veins, extra- ven- 
tricular CSF, lower arteries, lower capillaries, and 
lower veins, said non-vascular compartments including 
at least one of lower tissue, brain, ventricular CSF, and 
extra-ventricular CSF; 

establishing an atmosphere compartment, said atmo- 
sphere compartment representing a space located out- 
side the body; 

deriving a plurality of differential equations, each of said 
plurality of differential equations governing a pressure 
dynamic of one of said plurality of compartments and 
said atmosphere compartment; 
solving said plurality of differential equations; and 
producing a modeled circulatory system dynamic of the 
body based on said solving of said plurality of differ- 
ential equations. 

35. A mathematical model for simulating pressure dynam- 
ics of an intracranial system, comprising: 
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a means for dividing the body into a plurality of com- 
partments and a representation of a heart pump, each of 
said plurality of compartments representing a portion of 
the body, said plurality of compartments including at 
least one intracranial compartment and at least one 
extracranial compartment; 

a means for establishing an atmosphere compartment, 
said atmosphere compartment representing a space 
located outside the body; 

a means for deriving a plurality of differential equations, 
each of said plurality of differential equations go vein- 
ing a pressure dynamic of one of said plurality of 
compartments and said atmosphere compartment; 

a means for solving said plurality of differential equa- 
tions; and 

a means the outputting a modeled circulatory system 
dynamic of the body based on said solving of said 
plurality of differential equations. 
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